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ABSTRACT 

\ 

We  study  a  system  of  quasi linear  hyperbolic  conservation  laws  which  is 
hyperbolic  but  not  strictly  hyperbolic.  Such  systems  arise  naturally  in 
continuum  mechanics  such  as  elastic,  multiphase  flows.  We  are  interested 
mainly  in  the  large  time  behavior  of  the  solution.  Due  to  the  nonlinearity  of 
the  system  and  the  entropy  condition,  solutions  converge  to  very  simple 
elementary  waves.  Nonstrlct  hyperbolic! ty  of  the  system  may  cause  a  stronger 
nonlinear  interactions  between  waves  pertaining  to  different  families;  in 
particular,  such  interactions  may  regularize  linear  waves  in  the  solution. 

The  solutions  are  constructed  using  the  random  choice  method. 
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SIGNIFICANCE  AND  EXPLANATION 


We  study  a  system  of  quasilinear  hyperbolic  conservation  laws  which  is 
hyperbolic  but  not  strictly  hyperbolic.  Such  systems  arise  naturally  in 
continuum  mechanics  such  as  elastic,  multiphase  flows.  We  are  interested 
mainly  in  the  large  time  behavior  of  the  solution.  Due  to  the  nonlinearity  of 
the  system  and  the  entropy  condition,  solutions  converge  to  very  simple 
elementary  waves.  Nonstrict  hyperbollcity  of  the  system  may  cause  a  stronger 
nonlinear  interactions  between  waves  pertaining  to  different  families}  in 
particular,  such  interactions  may  regularize  linear  waves  in  the  solution. 

The  solutions  are  constructed  using  the  random  choice  method. 
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ON  A  HYPERBOLIC  SYSTEM  OF  CONSERVATION 
LAWS  WHICH  IS  NOT  STRICTLY  HYPERBOLIC 


Tai-Ping  Liu*  and  Ching-Hua  Wang** 


1.  INTRODUCTION. 


Consider  a  system  of  quasilinear  hyperbolic  conservation  laws 

n.D 

3t  3x 

where  U  -  U(x,t)  is  a  n-vector,  x  the  space  variable  and  t  the  time  variable.  One  of 
the  interesting  nonlinear  features  of  the  theory  of  conservation  laws  is  that  solutions 
attain  very  simple  asymptotic  states  as  t  tends  to  infinity.  The  system  is  strictly 
hyperbolic  if  3f(u)/3u  has  real  and  distinct  eigenvalues  X^u)  <  X2(u)  <  ...  <  Xn(u) 
for  each  state  u.  For  such  a  system,  the  solution  of  the  initial  value  problem  (1.1)  and 

U(x,0)  -  UQ(x>  (1.2) 

tends  to  elementary  waves  as  t  tends  to  infinity  [9]  ,  [10].  These  waveB  are  found  by 
solving  the  Riemann  problems  (1.1)  and 


U(x.O)  - 


u0(— )  for 

x  <  0 

U0(+*>)  for 

x  >  0 

(1.3) 


In  particular,  when  the  initial  data  (1.2)  have  a  compact  support  then  the  solution 
tends  to  the  zero  state,  [4]  ,  [2],  [9].  This  is  so  because  waves  combine  and  cancel  as  a 
consequence  of  the  nonlinearity  of  the  system  and  the  entropy  condition.  The  striking 
asymptotic  behavior  of  the  solution  can  be  understood  easily  for  scalar  conservation  law, 
[1],  [8],  [11]. 

When  the  system  is  non-strictly  hyperbolic,  that  is,  X^u)  may  equal  X^(u),  i  *  j, 
for  som  states  u,  then  waves  pertaining  to  different  characteristic  families  may  not 
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separate  as  tine  evolves*  When  this  happens,  nonlinear  interactions  of  these  waves  may 
alter  the  asymptotic  state.  In  this  paper  we  study  this  problem  for  a  system  of  two 
conservation  laws  whose  Riemann  problem  has  been  studied,  [5].  The  system  is  derived  from 
a  model  for  elastic  string.  One  of  the  characteristic  speeds  is  linearly  degenerate  and 
the  other  genuinely  nonlinear  in  the  sense  of  [6].  We  show  that  when  the  linear  wave  and 
the  nonlinear  wave  in  the  solution  of  (1.1)  and  (1.3)  are  separated  in  the  (x,t)-plane, 
then  the  asymptotic  behavior  of  the  solution  of  (1.1)  and  (1.2)  is  the  same  as  that  of  a 
strictly  hyperbolic  system.  In  this  case  the  asymptotic  state  consists  of  a  traveling  wave 
and  a  shock  or  rarefaction  wave,  c.f.  [10].  On  the  other  hand,  when  the  linear  wave  in  the 
solution  of  (1.1)  and  (1.3)  is  contained  in  the  nonlinear  wave  then  the  corresponding 
traveling  wave  becomes  substantially  more  regular  than  general  traveling  waves.  This  is  so 
because  of  the  strong  interaction  of  the  linear  and  nonlinear  waveB.  It  would  be 
interesting  to  investigate  the  problem  for  more  general  systems  where  the  interaction  of 
nonlinear  waves  of  different  families  occurs.  For  this  further  studies  are  necessary. 


2.  EQUATIONS  AND  RIEMANN  PROBLEM. 

The  following  two  conservation  laws  are  derived  from  a  model  of  elastic  strings,  [5] > 


(2.1) 


where  $  ■  $(u,v).  Let  (r,0)  be  the  polar  coordinates, 

2  2  2 

r  “  u  +  v  ,  tan  8  ■  v/u  , 
and  write  $  «  $(r,8).  The  characteristic  speeds  are 


h  -  * 


s  ■  ♦  * r  M  ■  i>?i 


(2.2) 


Thus  (2.1)  is  not  strictly  hyperbolic  on  J  3  {(r,8)  i  ■  0}  «  {(r,6)  i  1^  «  1^}.  The 
following  assumptions  on  $  are  consistent  with  physical  considerations. 
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<A)1  $(r,9)  ♦  »  as  r  ♦  0  or  r  * 

(A)2  *(r,8)  >  0, 

«>3  >  o,  i^>o, 


(A )  ♦(*,0)  la  convax  for  aach  fixed  9. 

4 

(A ) 1  and  (A)4  imply  that  ^  •  C,  C  any  constant,  la  a  simple  closed  curve.  Alao 
J  =  {(r,9,  -  -  0}  ia  a  simple  closed  curve.  (A)4  implies  that  >  X^  inside 

J  and  X 1  <  X^  outside  £.  The  first  characteristic  speed  X  ^  is  always  linearly 
degenerate  in  the  sense  of  [6].  (A)3  implies  that  X2  is  genuinely  nonlinear.  The  right 

eigenvectors  ,  i  -  1,2,  corresponding  to  X^  are  characterised  by 

V*  •  r,  -  0 
Vr  •  r2  -  0 

Along  the  curve  £,  where  X.)  ■  Xj,  the  system  is  diagonisable  when 
<A,5  If"*  0  on 

In  the  elastic  modal  it  is  reasonable  to  a  assume  that  $  is  a  function  of  r  only.  In 
this  «'ase  tne  above  assumptions  ( A)  1  ~  < A ) 5  are  satisfied  for  general  convex  $  with 
appropriate  growth  rate  at  r  ■  0  and  r  “  “. 

To  avoid  the  extreme  case  where  two  points  on  the  string  make  contact,  we  will  only 
deal  with  states  in  the  following  region i 

A  -  {0  i  r(0)  >  rQ,  1 6 ( u)  |  >i  0Q} 

for  a  fixed  rQ  >  0  and  0  <  0Q  <  ^.  The  Rlemann  problem  can  be  solved  by  similar  methods 
as  those  for  strictly  hyperbolic  systems.  We  have  three  kinds  of  elementary  waves.  A 
state  U  is  connected  to  U0  on  the  left  by  a  1-wave,  which  is  always  a  contact 


discontinuity,  if  U  c  T(uq)s 


T(Uq)  s  {0  i  *(U)  -  ♦iUQ)}  . 


(UQ,U)  forms  a  2-shock  wave  (2-rarefaction  wave)  if  U  e  Sj (UQ )  (U  c  (UQ  > ) » 


S2(UQ)  5  lo  j  9(0)  -  r<0)  <  r(nQ)} 

A2(U0)  5  {0  «  9(0)  -  8(Uq),  r(0)  >  r (UQ )} 


The  speed  of  a  2-shock  wave  (0a,0.)  is 


.  r1*1-r2*2 

o(u  ,u  )  -  -  ---- — 

12 


TO  solve  the  Riemann  problem  one  needs  to  know  how  the  speed  of  a  2-wave  Is  related  to  that 

of  a  1-wave.  Depending  on  the  relative  position  of  the  waves  in  the  solution  of  the 

Riemann  problem  (0^,0^)  we  have  the  following  cases: 

(I)  Ut  and  0^  are  both  outisde  J  and  $(0^)  >  $(0^). 

The  solution  consists  of  a  1-wave  (0,0  >  followed  by  a  2-rarefaction  wave 

Jt  w 

(0^,0^).  The  state  Um  is  characterized  by  $(On)  "  $(0^)  and  8(0^)  "  and  is 

outside  J. 

(IX)  Oj  and  are  both  inside  J  and  $(0^)  >  • 

The  solution  consists  of  a  2-rarefaetion  wave  (0,,0  )  followed  by  a  1-wave 

l  m 

<U  ,IJJ.  The  state  0_  is  inside  T  satisfying  8(0  )  -  8(U.)  and  <f>(U  )  -  $(0.,). 
n  A  ®  u  n  i  in 

(III)  Ut  is  inside  J  and  0^  is  outside 

The  solution  consists  of  a  2-rarefaction  wave  (0,,0  )  a  1-wave  (u  ,0  )  and  a 

*,  oi  m  n 

2-rarefaction  wave  (U  ,o_).  0  and  0  are  both  on  5"  satisfying  8(U  )  «  8(0.)  and 

njln  m  u  m  l 

8(0  )  -  8(0.,). 

The  above  three  cases  deal  with  solutions  containing  2-rarefactlon  waves.  The 

following  two  cases  deal  with  solutions  containing  a  2-shock  wave.  Thus  O^  is  assumed  to 

be  "closer"  to  the  origin  than  0^.  Let  0n  and  Un  be  defined  by: 

*<0  )  -  *(0),  8(U  )  -  8(0.)  , 
ni  /i  m  x 

*(0  )  -  $(U„>,  8(0  )  -  8(0.)  . 
n  1  n  K 

The  remaining  two  cases  are: 

(IV)  +(0^)  >  *(0^) 

The  solution  consists  of  a  1-wave  <0„,0  )  followed  by  a  2-shoek  wave  (0  ,0.).  It 

In  n  a. 

is  clear  that  o(0  ,0  )  >  4(0.)  which  is  the  speed  of  (0.,0  ). 

n  m  Jt  *  n 

(V)  <  *(0^) 

The  solution  consists  of  a  2-shock  wavs  (0.,0  )  followed  by  a  1-wave  (0  ,0*). 

k  m  in  't 
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3.  EXISTENCE  OF  THE  SOLUTION. 

The  Gil mm  scheme,  13],  for  strictly  hyperbolic  systems  can  also  be  used  to  construct 
solutions  for  general  hyperbolic  systems.  Choose  a  random  sequence  {a^)^,  -1  <  <  1  < 

and  mesh  length  Ax  ■  h,  At  “  s,  h/s  “  constant  satisfying  the  Courant-Friedrichs-Lewy 
condition: 

(C-F-L)  -  >  max  |X. (U) | 

8  u  1 

for  all  u  under  consideration.  The  approximate  solution  U^fx.t)  =  U.  (x,t»  {a.  })  is  a 
step  function  of  x  for  each  fixed  t  -  ks,  k  -  0,1,2,...,  with  possible  discontinuities 
at  x  “  jh,  j  +  k  ”  even.  By  resolving  these  discontinuities  (see  Section  2)  Ur(x,t)  is 
defined  for  t  e  (ks,(k  +  1)s).  Elementary  waves  issued  from  x  “  jr,  j  +  k  *  even,  do 
not  interact  before  t  -  (k  +  1)s  due  to  (C-F-L)  condition.  At 
t  «  (k+  1)s,  Ur(x,t)  is  not  a  step  function  and  the  random  sequence  is  u8e<J 

to  approximate  it  by  a  step  function: 

Uh(x1(k+1)s  +  0)  s  Uh((j  +  ak)hft(k+1)s  -  0),  (j-1)h  <  x  <  (j+1)h,  j  +  k  -  odd  . 

This  defines  inductively  the  approximate  solution  Uh!x,t)  for  all  t. 

The  convergence  of  the  approximate  solution  (x, t ) }  as  h  tends  to  zero,  is 

proved  in  two  steps:  First,  one  shows  that  U^fXjt)  has  bounded  variation  in  x  for  each 

fixed  t.  This  implies  by  diagonal  process  that  {Uh(x,t)J  converges  strongly  in  L^x) 

for  any  rational  t.  To  prove  the  convergence  for  all  t,  one  needs  the  Lipschitz 

continuity  of  L.,(x)  in  t.  This  follows  from  the  estimate  on  the  total  variation  in  x 

of  U.  and  that  U.  (x,t)  has  a  finite  speed  of  propagation.  For  this  the  system  has  only 
h  n 

to  be  hyperbolic,  not  necessarily  strictly  hyperbolic.  Thus  for  the  existence  of  the 

solution  one  need  only  to  estimate  the  total  variation  in  x  of  the  approximate  solution 

U.  (x,t).  It  is  convenient  to  introduce  a  new  coordinate  (r,6)  as  follows:  For  any  U 
n 

in  the  region  A  let  U  be  the  unique  state  on  the  same  side  of  J  with  $(U)  *  $(U) 
and  0(U)  “  0.  We  set 

r(U)  «  r (U)  . 
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It  follows  from  the  assumptions  ( A ) 1  ~  ( A)  ^  of  Section  2  that  the  transformation 

(r,0)  -  (u,v)  Is  nonsingular  for  U  -  (u,v)  in  a  bounded  region  of  A.  Note  that  a  1- 

wave  takes  values  along  r  -  constant  and  a  2-wave  takes  values  along  6  “  constant.  The 

strength  of  a  1-wave  is  defined  by  the  jump  of  8  across  the  wave;  and  the  strength  of  a 

2-wave  by  the  jump  of  r  across  it.  The  following  lemmas  follow  directly  from  the  recipe 

for  solving  the  Riemann  problem  presented  in  the  last  section.  We  omit  details. 

Lemma  1 ;  Suppose  that  the  initial  data  U( x, 0 )  stay  in  a  bounded  region 

fl  =  {u  i  0  <  r1  <  r(U)  <  r2  <  »,  |8|!u)  <  8Q  <  |}. 

Then  any  approximate  solution  also  stay  in  ft  for  any  (x,t). 

Lemma  2:  Suppose  that  U_ ,U  ,U „  are  three  states  in  ft,  then  the  total  strength  of  waves 
■  In'1 

in  the  Riemann  problem  (17  ,U^)  is  no  larger  than  the  sum  of  the  total  strength  of  waves 
in  the  Riemann  problem  and  the  total  strength  of  waves  in  the  Riemann  problem 

<w- 

Lemma  1  shows  that  there  are  bounded  invariant  regions  for  the  solution  and  so  in 
particular  the  (C  -  P  -  L)  can  easily  be  satisfied.  Lemma  2  shows  that  nonlinear 
interactions  do  not  cause  an  increase  in  the  strength  of  waves.  This  is  sufficient  to 
obtain  the  desired  estimate  on  the  variation  in  x  of  the  approximate  solution  U^tx.t), 
(c.f.  [13]).  Thus  we  have  the  following  existence  theorem. 

Theorem:  Suppose  that  assumptions  (A) ^  ~  (A)g  hold  and  the  initial  data  U(x,0)  stay  in 
a  bounded  region  ft  =  (u  :  0  <  r1  <  r(U)  <  r2  <  »,  |8(U)|  <  e0  <  f}  and  have  bounded 
variation  locally  in  x.  Then  the  initial  value  problem  for  (2.1)  has  a  global  solution 
U(x,t)  which  stays  in  ft  and  has  bounded  variation  locally  in  x. 

4.  ASYMPTOTIC  BEHAVIOR  OF  THE  SOLUTION. 

In  this  section  we  assume  that  w  “  ^(r).  The  advantage  of  this  assumption  is  the 
weaker  coupling  of  the  system  (2.1)  and  so  the  behavior  of  2-waves  can  be  studied 
Independently.  Across  a  1-wave,  the  value  of  r  is  unchanged.  A  2-rarefactlon  wave  has 
speed  (r$)r  and  a  2-shock  wave  (Ug,^)  has  speed 
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.  riJ|i  ro*o 

Moreover,  across  a  2-wave  9  is  constant.  Therefore  the  behavior  of  r,  and  hence  the 
behavior  of  2-waves,  can  be  described  by  the  scalar  conservation  law 

||  .  JUJlEii  .  .  . 

A  1-wave  propagates  with  speed  $  and  so  1 -waves  are  described  by 

i  +  »(r>  -M  -  0  •  <4-2> 

The  identities  (4.1)  and  (4.2)  can  be  derived  from  the  system  (2.1)  directly.  For 
instance,  multiplying  the  first  equation  of  (2.1)  by  u  and  the  second  equation  by  v  and 
summing  them  up,  one  obtains  (4.1).  This  procedure  is  justified  for  smooth  solutions 
(u,v).  However,  in  general  when  an  equation  is  derived  from  a  system  of  conservation  laws 
through  nonlinear  transformations,  the  derived  equation  may  not  be  a  consequence  of  these 
conservation  laws  in  the  weak  sense.  In  the  present  situation  the  procedure  is  justified 
because  the  jump  condition  and  entropy  condition  for  (4.1)  and  (4.2)  are  consistent  with 
those  for  (2.1).  Our  strategy  is  to  study  the  behavior  of  2-waves  UBing  (4.1)  and,  having 
obtained  the  behavior  of  r,  use  (4.2)  to  study  the  behavior  of  1-waves. 

From  now  on  we  assume  that  the  initial  data  satisfy 


U(x,0) 


U.  for  x  <  -s/2 
t 

U„  for  x  >  s/2  . 
A 


for  some  constant  states  U.  and  U,.  We  denote  by  U(x,t)  the  solution  of  the 

l  A 

corresponding  Riemann  problem  (2.1)  with 


U(x,0) 


for  x  <  0  , 
0^  for  x  >  0  . 


The  behavior  of  2-waves  is  described  by  (4.1)  with 


r(x,0) 


r.  for  x  <  -s/2 
l 

r„  for  x  >  s/2  , 

A. 


There  two  cases;  r 8 ) : 
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Case  Is  r.  >  r„ 

Jt  A 

There  exists  T  >  0  such  that  for  t  >  T,  the  solution  r(x,t)  of  (4. 

is  a  shock  wave  (r.,r): 

Z.  A 

r  for  x  -  x  <  at 

r(x,t)  -  *  0  ,  t  >  T  . 

r.  for  x  -  x„  >  at 
A.  0 


V(V  -  V(ri> 
^ — 


Xn  =  - 

0  r.  -  r 


9/2  r  +  r 

/  [ - - - r(x,0)  ]dx  . 


1  -s/2 


Case  2:  r  < 

Z  A 

The  solution  r(x,t)  tends  to  a  generalized  N-wave  defined  as  follows: 


N(x,t)  =  N(x,t;p,q;r 4,r t) 


f  for  x  <  XjCr^Jt  -  /^pA'Cr^t  £  x(t) 


for  x  >  Xj(r ^)t  +  /2qX^(r^)t  =  x^tt) 
r(C)  otherwise:  5  -  x/t 


p  =  min  /  [r ( y ,0 )  -  r  Jdy,  q  =  max  f  [r(y,0)  -  r  1  dy 

X  -»  *  XX 


More  precisely,  there  exist  Lipschitz  continuous  curves  x  *  x^(t)  and  x^(t 
(-s/2,0)  and  (s/2,0),  respectively,  with  the  following  properties: 

(i)  Ix^t)  -  x^(t)|  +  !xi(t)  -  x  (t>|  -  0(S). 

(il)  |r(x,t)  -  N(x,t)|  “  0( 1 )t  1  for  (x,t)  between  x„(t)  and  x„( 

l  A 

between  x.(t)  and  x„(t). 

I  A 

(iii)  |r(x,t)  -  N(x,t)|  »  0(1)t  either  <x,t)  lies  between  x^ft) 
or  between  x^ft)  and  x^(t). 

( iv)  r(x,t)  »  r^  for  x  <  x^(t)  , 

r^  for  x  >  x^(t)  . 


1)  and  (4.3j’ 


)  through 


)  and  also 


and  x^(t) 
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We  now  use  the  above  known  results  to  study  the  behavior  of  1 -waves.  When  waves  of 
different  families  in  the  solution  of  the  Riemann  problem  (2.1)  and  (4.3)  are  separated, 
the  asymptotic  behavior  of  the  solution  U(x,t)  of  (2.1)  and  (2.2)  is  similar  to  that  of 
solutions  of  a  strictly  hyperbolic  system.  In  this  case  1 -waves  tend  to  a  traveling 
wave.  We  illustrate  this  by  investigating  Case  (V)  of  Section  2.  Thus  we  have 
<  4>(r^).  From  (4.2)  characteristic  curves  for  1 -waves  are  given  by: 

“  =  $(r(x,t ) )  . 

Through  (-s/2,0)  and  (s/2,0),  respectively,  we  draw  two  characteristics  and  X2» 

We  know  from  the  above  discussion  of  2 -waves,  (case  1),  that  r(x,t)  is  a  shock  wave  for 
t  >  T.  The  speed  of  the  shock  wave  is  a  which  is  less  than  4>(r^)  because 
4>(r  )  <  4>(r^).  Consequently  the  1-characteristic  X1  lies  to  the  right  of  the  shock  wave 
after  time  T^,  T^  finite.  After  time  T1  there  exists  no  1-wave  either  to  the  left 
of  X ^  or  to  the  right  of  X2,  and,  between  and  X2,  0  is  constant  along 

dx/dt  *  ^(r^)  according  to  (4.2).  Thus  after  time  Tj,  1 -waves  become  a  traveling  wave 
between  and  X2  and  with  speed  $(r^).  The  value  of  S  to  the  left  of  X1  is  9^ 

and  to  the  right  of  X2  is  0^. 

The  situation  is  more  complicated  when  waves  of  different  families  do  not  separate. 

We  exemplify  this  by  investigating  Case  III  in  Section  2.  In  this  case,  the  solution  of 


the  Riemann  problem  (U  .u)  is  a  contact  discontinuity  (U  ,U  )  sandwiched  by  two 
Jc  a  m  n 

2-rarefaction  waves  (U„,U  )  and  (U  ,U.).  The  states  Um  and  are  on  V  and  so 

£  m  n  h.  m  n 

A  „  (U  )  =  A.(U  )  =  A„(U  )  =  A.(U  ).  We  consider  the  case  when  U  t  U  and  U  *  UH. 

1m  2  m  in  2  m  £  m  n  a 

Note  that  A^  attains  an  absolute  minimum  on  T  and  so  the  characteristic  curve  for  (4.2) 

always  has  speed  larger  than  or  equal  to  ^(r  )  =  <{>(r  ).  Suppose  that 
x 

e(x)  =  f  [r  (y)  -  r  )dy  attains  a  minimum  at  x  =  xn.  Then  the  characteristic  line 
i  0  m  u 

0 

through  (xQ,0)  for  (4.1)  exists  for  all  t  >  0,  (7).  In  this  case  the  ^naracteristic 

curve  through  (xQ,0)  for  (4.2)  coincides  with  that  for  (4.1).  In  fact  only  such 

characteristic  curves  for  (4.2)  may  travel  with  the  minimal  speed  $(r  )  for  all  t  >  0. 

m 

Thus  when  I(x)  takes  minima  at  x  -  a'  and  x  =  a"  and  not  for  any  x  e  (a', a"),  then 
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* 


all  characteristic  curves  through  (x,0),  a'  <  x  <  a"«  tend  to  the  characteristic  line 


x  =  x.  +  4(r  )t  and  so 
0  Y  m 


lim  9<x  +  $(r  )t,t)  »  9(a'  +  0,0)  for  x  e  (a',aM) 
^  m 


(4.4) 


It  is  clear  that 


m(r  )  =  (x.  :  /  [r„(xl  -  r  Jdx  attains  minimum  at  x  «  x. } 

m  0  o  0  m  ” 

is  a  closed  set.  Denote  by  xffl  and  xM,  respectively,  the  smallest  and  largest  numbers 

in  m(r  ).  The  set  Mr  )  =  [x  ,xl  -  m(r  )  is  an  open  set.  We  have  just  described  the 

ni  m  m  M  m 

asymptotic  behavior  of  9(x,t)  for  x  =  x  +  *(r  )t,  x  in  a  component  of  Mr  ).  When 

u  in  u  m 

m(r  )  contains  an  interval  (a,b)  we  have 


r(x,t)  =  r.  for  *(r  )t  +  a  <  x  <  *(r  )t  +  b  , 

Z  ra  m 

0(x,t)  =  0(x  -  4(r  )t,0),  for  4(r  )t  +  a  <  x  <  Mr  )t  +  b 
m  m  m 


(4.5) 


For  x  <  x  +  Mr  It,  all  characteristic  curves  for  (4.2)  tend  to  the  characteristic  line 
m  m 


through  (x  ,0).  Since  0(x,O)  -  9„  for  x  <  -3/2  we  have 

m  z 

lim  0(x  +  Mr  )t,t)  »  0.  for  x  <  x  . 

ra  l  m 


(4.6) 


We  have  seen  that  the  initial  values  of  9  restricted  to  (-«°,x„)  may  not  be  carried  to 
t  =*  «°  and  the  asymptotic  shape  of  9(x,t),  x  -  4(r  )t  e  (-<»,r  ),  is  in  general  a  step 
function.  On  the  other  hand,  the  initial  data  b(x,0),  x  >  M,  are  in  general  carried  to 
t  =  «■  and  the  asymptotic  behavior  of  0(x,t),  x  -  Mr  )t  e  (x  „,<»),  is  a  traveling  wave 
taking  values  (0(x,O),  x  >  x  > .  We  now  show  that  this  traveling  wave  has  a  finite 
width.  For  this,  we  draw  a  1-characteristic  X  through  (s/2,0)  and  estimate  the 
asymptotic  location  and  speed  of  X.  Since  r(x,t)  tends  to  an  N-wave  and  all  2-shock 
waves,  except  the  one  issued  from  s/2,  decay  at  the  rate  1/t,  we  will  carry  out  the 
analysis  by  supposing  that  r(x,t)  is  a  centered  rarefaction  wave.  This  is  done  for 
simplicity;  the  general  case  can  be  dealt  with  similarly.  Thus  it  follows  from  the 
structure  of  centered  rarefaction  wave  that  when  r  =  r^  at  {t  «  t^}  n  X  and 
r  =  -  Ar  at  (t  ■*  t^  +  At}  n  X  then 

U2(r)  -  X2<r  -  Ar))t  «  U2(r)  -  X^rllAt 
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This  yields  an  ordinary  differential  equation  for  the  value  of  r  along  X: 


dl2(r(t)) 

dt 


■U2<r<t))  -  X^rttni/t 


(4.7) 


Note  that  X2(r)  >  X^r)  for  r  in  the  region  under  consideration,  i.e.  r  >  rB.  Thus 


the  above  identity  implies  that 


r(t)  ■*  r  as  t  + 
m 


dg(t) 

dt 


The  rate  of  this  convergence  is  determined  as  follows:  Consider 

-C<t)/t 

5<t>  =  Xj  (r(t) )  -  X^rtt))  . 

It  is  clear  that 


?(t)  »  0(  1)  £  as  t  ♦ 


Note  that  as  r(t)  tends  to  r_ 


I  X, ( r ( t )  )  -  X,(r  )|  ~  |r(t)  -  r  I 
z  z  m  in 


I  X  (r( t) )  -  X,(r  ) I  -  |r(t)  -  rl 
1  l  m  m 


(4.8) 


and  so  r(t)  and  £(t)  have  the  same  qualitative  behavior.  We  thus  have: 


!r(t)  -  r  |  »  0(1)  ±  , 
in  t 


IX. (t)  -  X, ( r  >|  -  0(1)  —  . 

i  i  n\  ^z 


(4.9) 


This  implies  that  X  tends  to  a  straight  line  x  -  4>(r  )t  “  constant  at  the  rate  1/t.  In 

m 

particular  0(x,t),  x  >  ^  +  4> ( ) t ,  tends  to  a  traveling  wave  of  finite  width.  Note  that 

the  speed  of  X  is  always  larger  than  4(r  )  and  so  the  distance  between 

in 

{(x,t)  :  x  «  x„  +  *(r  )t}  and  X  is  an  increasing  function  of  time.  Consequently, 

M  m 

except  for  the  exceptional  case  where  XM  m  S/2,  the  asymptotic  distance  of  these  two 
curves  is  finite  and  positive.  This  completes  the  description  of  the  asymptotic  state  of 
9(x,t) . 
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In  the  above  arguments  we  have  assumed  that  8(x,t)  is  constant  along  the 
2-characteristics  dx/dt  «  $(r(x,t))  and  that  2-characteristics  are  defined  for  all 
t  >  0.  This  is  so  because  across  a  2-shock  wave  6  is  unchanged  and  the  1-characteristic 
speeds  $  on  both  sides  of  the  shock  are  either  greater  or  less  than  the  shock  speed. 

The  above  arguments  can  be  applied  to  treat  other  cases  in  section  2.  We  briefly 
state  the  asymptotic  results  in  the  following  theorem. 

Theorem  4.1;  Suppose  that  is  a  function  of  r  and  hypothesis  (A)^,  i  -  1,2,... ,5, 
hold  and  the  initial  data  U(x,0)  equals  for  x  <  -S/2  and  for  x  >  S/2.  Then 

the  asymptotic  behavior  of  the  solution  of  (2.1)  is  as  follows: 

(I)  The  behavior  of  2-waves  is  described  by  (4.1).  Thus  2-waves  tend  to  a  single  2- 

shock  wave  when  r,  >  and  to  a  N-wave  when  r  <  r„ . 

I  A.  I  A 

(II)  When  1 -waves  and  2-waves  in  the  solution  of  the  Riemann  problem  (U^ ,11^)  for 
(2.1)  are  separated,  1-waves  tend  to  a  traveling  wave  of  finite  width  in  finite  time.  The 
traveling  wave  assumes  all  the  values  of  0(x,O).  The  same  also  holds  when  the  2-wave  in 
the  solution  of  (U^ ,U^)  is  a  shock  wave. 

(iii)  When  the  2-wave  in  the  solution  of  is  8  rarefaction  wave  which 

contacts  the  1-wave  in  ( U ^ ) ,  1-waves  tend  to  a  traveling  wave  as  t  tends  to 
infinity.  The  traveling  wave  is  a  combination  of  a  step  function  and  a  general  traveling 
wave.  Moreover,  the  traveling  wave  does  not  assume  all  the  initial  values  0(x,O)  and  is 
described  by  (4.4),  (4.5)  and  (4.6). 
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